NASA Contractor Report 4763 


Evaluation of the MV (CAPON) Coherent 
Doppler Lidar Velocity Estimator 

B. Lottman and R. Frehlich 

Cooperative Institute for Research in Environmental Sciences (CIRES) 

University of Colorado • Boulder, Colorado 



National Aeronautics and Space Administration Prepared for Marshall Space Flight Center 

Marshall Space Flight Center • MSFC, Alabama 35812 under Grant NAG8-253 


January 1997 




Contents 


1 Introduction 1 

2 Typical Parameter Regimes 1 

3 Velocity Estimator Performance 3 

4 Potential Uses 4 

5 Computational Effort 4 

6 Results M=10 5 

7 Results M=16 9 

8 Results M=32 10 

9 Results M=64 12 

10 Results M=128 13 

11 Results M=512 14 

12 Acknowledgements 16 

Appendix A, M=16 17 

Appendix B, M=32 21 

Appendix C, M=64 27 

Appendix D, M=128 32 

Appendix E, M=512 37 


iii 



List of Tables 


Table 1. Parameter regime for a typical 2pm and 10 pm LAWS type lidar. 2 

Table 2. Parameter regime for a typical 2/im and \0pm boundary layer lidar. 2 

Table 3. Parameter regime for typical 2/im and 10 pm LAWS type lidar observing a 
thick cloud. 3 

Table 4. Parameter regimes evaluated in this paper. 3 

Table 5. Computational effort for the CAPON Estimator based on estimator order 
p, oversampling os, and the number of data points in the range cell M . 4 

Table 6. Model Fit Parameters for Af=10, f7=0.2 and 0.3. 8 

Table 7. Model Fit Parameters for M= 16, $7=0.2, 0.3, 0.5 and 1.0. 10 

Table 8. Model Fit Parameters for M— 32, £7=0.2, 0.3, 0.5, 1.0, 2.0, and 3.0. 11 

Table 9. Model Fit Parameters for M=64, $7=0.3, 0.5, 1.0, 2.0, and 3.0. 13 

Table 10. Model Fit Parameters for M= 128, $7=0.5, 1.0, 2.0, 3.0, and 7.0. 14 

Table 11. Model Fit Parameters for M=512, $7=7.0 and 10.0. 15 

List of Figures 

Figure 1. Estimator performance as a function of oversampling os and order p for 
M=32 and $7=0.5. 3 

Figure 2. Estimator performance vs. computational effort for M= 512 and $7=7.0. 5 

Figure 3. Performance vs. oversampling os for M=10 and $7=0.2. 6 

Figure 4. Performance vs. oversampling os for M=10 and $7=0.3. 6 

Figure 5. Performance vs. order p for the CAPON estimator for M— 10 and $7=0.2. 6 

Figure 6. Performance vs. order p for the CAPON Estimator for M=10 and $7=0.3. 6 

Figure 7. Single realizations of the CAPON estimates for M — 10 and $7=0.2 for 
$=1,10,50,1000 with optimal estimator parameters. 7 

Figure 8. Single realizations of the CAPON estimates for M=10 and $7=0.3 for 
$=1,10,50,1000 with optimal estimator parameters. 7 

A 

Figure 9. PDF of the Velocity Estimates using the CAPON estimator with M = 10 
and $7=0.2 for $=1,10,50,1000. 7 

iv 



Figure 10. PDF of the Velocity Estimates using the CAPON estimator with M=10 
and £2=0.3 for for $=1,10,50,1000. 7 

Figure 11. Performance of the CAPON estimator with M=10 and £2=0.2 and 0.3 
using optimal estimator parameters. Best fit model is (-). 8 

Figure 12. Optimal parameters for the CAPON estimator with M=10 and £2=0.2 
and 0.3. g 

Figure 13. Minimum estimator order p with M = 10 and £2= 0.2 and 0.3 for no 
greater than 5the good estimates. g 

Figure 14. Performance of the CAPON estimator with M — 16 and £2=0.2, 0.3, 0.5, 
and 1.0 using optimal estimator parameters. 9 

Figure 15. Optimal parameters for the CAPON estimator with A/=16 and f2=0.2, 

0.3, 0.5, and 1.0. 9 

Figure 16. Minimum estimator order p with M= 16 and £2=0.2, 0.3, 0.5, and 1.0 for 
no greater than 5% and 10% increase in the standard deviation g of the good 
estimates. q 


Figure 17. Performance of the CAPON estimator with M = 32 and £2=0.2, 0.3, 0.5, 

1.0, 2.0, and 3.0 using optimal estimator parameters. 


10 


Figure 18. Optimal parameters for the CAPON estimator with A/=32 and £2=0.2, 
0.3, 0.5, 1.0, 2.0, and 3.0. H 

Figure 19. Minimum estimator order p with Af= 32 and £2=0.2, 0.3, 0.5, 1.0, 2.0 
and 3.0 for no greater than 5% and 10% increase in the standard deviation g of 
the good estimates. jj 

Figure 20. Performance of the CAPON estimator with M=64 and £2=0.3, 0.5, 1.0, 

2.0, and 3.0 using optimal estimator parameters. 12 

Figure 21. Optimal parameters for the CAPON estimator with M=64 and £2=0.3, 
0.5, 1.0, 2.0, and 3.0. 12 

Figure 22. Minimum estimator order p with M — 64 and £2=0.3, 0.5, 1.0, 2.0 and 
3.0 for no greater than 5% and 10% increase in the standard deviation g of the 
good estimates. 12 

Figure 23. Performance of the CAPON estimator with M=128 and £2=0.5, 1.0, 2.0, 

3.0, and 7.0 using optimal estimator parameters. 13 

Figure 24. Optimal parameters for the CAPON estimator with M=128 and £2=0.5, 

1.0, 2.0, 3.0, and 7.0. 14 


v 



Figure 25. Minimum estimator order p with A/=128 and 0=0.5, 1.0, 2.0, 3.0 and 
7.0 for no greater than 5% and 10% increase in the standard deviation g of the 
good estimates. 

Figure 26. Performance of the CAPON estimator with M=512 and $1=7.0 and 10.0 
using optimal estimator parameters. 

Figure 27. Optimal parameters for the CAPON estimator with M=512 and $1=7.0 
and 10.0. 

Figure 28. Minimum estimator order p with M=512 and $1=7.0 and 10.0 for no 
greater than 5% and 10 % increase in the standard deviation g of the good 
estimates. 


vi 



1 Introduction 

The first coherent Doppler lidars were based on 
the CO 2 laser [1, 2. 3, 4]. More recently, solid 
state Doppler lidars have been successfully op- 
erated [5, 6, 7]. Coherent Doppler lidar is un- 
der consideration for measurements of the global 
wind field from space [[3, 8, 9]]. A coherent 
Doppler lidar operating at wavelength A can es- 
timate the radial component of the velocity v = 
A//2 as a function of range for every transmitted 
pulse by estimating the Doppler frequency shift 
/. The Probability Density Function (PDF) pro- 
vides a complete statistical description of an es- 
timator. For the better velocity estimators, the 
PDF is characterized by a fraction b of uniformly 
distributed bad estimates or random outliers and 
a distribution of good estimates centered on the 
true mean velocity. Frehlich and Yadlowsky[10] 
investigated the performance of mean-frequency 
estimators by modeling the distribution of good 
estimates as a Gaussian PDF with standard de- 
viation g . Approximately universal curves of 
performance of velocity estimates (and mean- 
frequency estimates) are produced by plotting 
the standard deviation g normalized by the sig- 
nal spectral width in velocity space w v — Xw/2 
versus the parameter $ = SNR • M for fixed 
ft = wMT$ where w is the signal spectral width, 
Ts is the sampling interval of the complex data, 
SNR is the Signal-to-Noise Ratio, and M is the 
number of complex data points per observation. 
ft is fixed for fixed range resolution and fixed 
transmitted pulse length, the typical comparison 
case for lidar performance. For shot-noise domi- 
nated operation, the parameter <J> = t}hN v where 
iff! is the heterodyne efficiency (77// ~ 0.4 in the 
far field. The far field condition is usually valid 
for space- based platforms.) and N e is the aver- 
age number of photo-electrons per estimate. <f> is 
the average number of coherent photo-electrons 
per estimate. The parameter ft is proportional 
to the number of independent samples of the sig- 
nal per range gate, sometimes called the “speckle 
count,’'. 

Anderson [11] proposed the Minimum Vari- 
ance (MV) or CAPON spectral estimator for 


space-based measurements of winds. Assuming 
and ft are known, we evaluate the perfor- 
mance of the CAPON estimator where the auto- 
regressive coefficients up to order p are deter- 
mined from the biased covariance estimates and 
the Yule- Walker solution [10. 12, 13]. The op- 
timal order p and the optimal sampling interval 
for estimating the maximum of the spectral esti- 
mate S(f) for the CAPON estimator are deter- 
mined by minimizing the standard deviation of 
the good estimates. Various parameter regimes 
will be considered. 

2 Typical Parameter Regimes 

A variety of parameter regimes relevant to space- 
based measurements are evaluated in this paper. 
The basic space-based parameters were deter- 
mined by the following criteria; 

1) A velocity search space v s of 50 m/s was se- 
lected. The sampling interval T s for complex 
data is given by T s = X/(2v s ). The sampling 
interval for 2 and 10 pm lidar is 20 and 100 ns 
respectively. 

2) For numerical weather prediction, current 
models use a 1 km height resolution. This im- 
plies a 1.5 km range cell (line of sight) for typical 
satellite scanning geometries. The range resolu- 
tion is Ap-h Ar where Ap is the distance the pulse 
moves per velocity estimate, i.e., Ap = MT s cj 2 
where c [m/s] is the speed of light. Ar [m] is 
the full width half max (FWHM) spatial extent 
of the pulse, i.e. Ar = Atc/2 where A t [s] is 
the FWHM temporal extent of the pulse. A A;; 
of 1500m was selected. Since M = 2Ap/(T s c), 
M — 500 and 100 for 2 and 10 pm lidars re- 
spectively. For a large Ap the spectral width 
w of the returned signal is usually dominated 
by wind turbulence. The lidar signal is approx- 
imated as a Gaussian with a spectral width in 
velocity space w v of 1 m/s and since w = 2w v /X, 
w is 1 and 0.2 MHz for a 2 and 10 pm lidar re- 
spectively. ft = wMT s where ft is proportional 
to the speckle count, ft = 10.0 and 2.0 for a 2 
and 10 pm lidar respectively. 

Examples of parameter regimes for a 2pm and 
1 Op, 77i LAWS type lidar are shown in Table 1. 
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A p 

M 

n 
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20 

1.0 

1536 

512 

10.2 

2 

20 

1.0 

768 

256 

5.1 

2 

20 

1.0 

384 

128 

2.5 

10 

100 

0.2 

1920 

128 

2.5 

10 

100 

0.2 

960 

64 

1.3 

10 

100 

0.2 

480 

32 

0.6 


Table 1. Parameter regime for a typical 2 /xm 
and 10 fim LAWS type lidar. A is in /xm . T s is in 
ns. w is in MHz. A p is in m. 


Typical lidar parameters for boundary layer 
measurements are also of interest. The basic pa- 
rameters for the typical boundary layer lidar are 
based on the following; 

1) 50 m/s was selected for the velocity search 
space. 2 and 10 /xm was selected for lidar wave- 
length A. The sampling time T s for complex data 
for a 2 and 10 /xm lidar is 20 and 100 ns respec- 
tively. 

2) The boundary layer is typically 1000 m in 
depth. To characterize activity within the layer, 
measurements are usually taken at least every 
100 m. A range cell length A p of 48m was se- 
lected and M = or M = 16. 

-L S C 

For a small A p the spectral width w of the 
returned signal is determined by the transmit- 
ted pulse (pulse dominated) and not wind tur- 
bulence. For a 2 /xm example, the FWHM 
pulse width A* is 0.3 /xs . The spectral width 
w (pulse dominated) is w = \/ln 2/2/7 tA t or 

- 0.1873/A* and ft = wMT s or ft = 0.2. 

For very large spectral width w , the velocity 
search space becomes comparible with w. In this 
regime the PDF’s do not have a uniform layer 
of outliers and the two parameter model is not 
valid. The upper bound on ft is when the 16a 
spectral width of the signal approachs the full 
velocity search space. For example, for M = 
16, 16ic < 1/T S , or wT s < 1/16, therefore ft = 
wMT, < M/16, or ft < 1. 

Examples of typical parameter regimes for the 
2/ini and 10/im boundary layer lidar are shown 
in Table 2. 

For a LAWS system, observation of thick 
clouds is of interest. Thick clouds generate a sig- 
nal similar to a hard target return and the signal 


spectral width is determined by the transmitted 
pulse. 

The basic parameter set for thick cloud obser- 
vation using a LAWS type system is determined 
as follows; 

1) Sampling interval is selected based on a veloc- 
ity search space of 50 m/s and the lidar wave- 
length A. The sampling interval T s for 2 and 10 
/xm lidar is 20 and 100 ns respectively. 

2) The observation time MT S is determined as a 
fraction of the FWHM of the transmitted pulse 
At. MT S = kAt ; where k is a user defined con- 
stant, which is related to ft by ft = ky/\n2//2 / tt 
or ft=.1873 k. At determines M. 
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Table 2. Parameter regime for a typical 2/xm 
and 10/im boundary layer lidar. A is in m. T s is 
in ns. At is in fis. A p is in m. 


This regime is similar to the high SNR small 
range gate scenario for boundary layer measure- 
ments. The parameter sets for the thick cloud 
target with various transmitted pulse widths 
for 2/xm (T s = 20ns) and 10/im (T 5 = 100ns) ar(? 
shown in Table 3. 

Tables 1 through 3 give the general regimes of 


2 





interest. The parameter sets evaluated in this 
paper are shown in Table 4. 
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Table 3. Parameter regime for typical 2 pm and 
10 pm LAWS type lidar observing a thick cloud. 
At is in [is. 



Table 4. Parameter regimes evaluated in this 
paper. 


3 Velocity Estimator Perfor- 
mance 

Simulations of complex coherent Doppler lidar 
data were produced by the method described 
in Frehlich and Yadlowsky[10]. This produces 
data with a specified auto- covariance function 
which is chosen as a Gaussian function defined 
by the parameters <J>, ft, and M [see Eq, (22) of 
Ref. [10]]. This is a good approximation for 2 pm 
coherent lidars[14]. 

For the CAPON estimator, performance is a 
function of both estimator order and oversam- 
pling. Oversampling refers to sampling the spec- 
tral function at smaller intervals to improve; es- 


timator accuracy. In velocity space, the sam- 
pling interval A?; = v s /(M os) where os is the 
oversampling, os = 1 means that the spectral 
function is sampled with M data points. The 
maximum value of the spectral estimate and the 
two closest neighbors are fit to a parabola. The 
velocity that produces the peak of the parabola 
is selected as the estimate. The value of x that 
maximizes f(x) for uniformly sampled x* is 

X t -Ti+ /( J 2) - /(g o) 

x peak~ *1 + 2{ 2 f {xi) -f( X0 )-f(x 2 )) 

where xq is the index of the spectral function to 
the left of the peak, x\ is the index of the spec- 
tral function at the peak, and X 2 is the index of 
the spectral function to the right of the peak. 

Generally, performance is less sensitive to 
oversampling than estimator order. Examples of 
oversampling and order sensitivity for M = 32 
and ft = 0.5 are shown in Fig. 1. Performance is 
normalized by the minimum standard deviation 
of the good estimates for each regime g mtn . 



Figure 1 . Estimator performance as a function 
of oversampling os and order p for A/=32 and 
SI— 0.5. Curves are offset by 0.10. The la error 
bars are less than the symbol size. 

The technique for determining the optimal es- 
timator parameters (p, os) for given basic pa- 






rameters (M . ft, $) is to first determine opti- 
mal oversampling os with an approximate opti- 
mal order, then determine optimal order p us- 
ing the oversampling result os. Parameters that 
give the best performance within the statistical 
limitations of the simulation are considered the 
optimal parameters. 

For large M with large ft, estimator accuracy 
does not improve with oversampling {os = 1 pro- 
duces the best performance). For the M = 512 
ft = 7.0 and 10.0 cases order sensitivity is com- 
puted but oversampling sensitivity is not com- 
puted. 

The results of each parameter set include: per- 
formance as a function of oversampling and or- 
der, single realizations of the CAPON spectral 
estimates, and PDF’s for the velocity estimators 
for a subset of cases using optimal estimator pa- 
rameters. 


4 Potential Uses 

The standard deviation of the good estimates 
g and the fraction b of uniformly distributed 
bad estimates can be used to generate simu- 
lated velocity data for a given parameter regime. 
To generate velocity data, velocity values are 
drawn from one of two random number gener- 
ators. The good estimates are represented by a 
Gaussian random variable with standard devia- 
tion g , while the bad estimates are represented 
by a uniform random variable over the velocity 
search space v s . Selection of a good or bad es- 
timate is based on the parameter b, the fraction 
of uniformly distributed bad estimates. 

For a LAWS type mission, the final data prod- 
uct will be generated from multiple shots in a 
resolution cell with a specified scanning configu- 
ration e.g. a conical scan[9] . The desired data 
product is the average wind vector over the res- 
olution cell in 3D. The PDF of any multishot 
velocity algorithm can be derived from the sin- 
gle shot PDF’s using a Monte Carlo analysis. 


5 Computational Effort 


Computational effort for complex data is shown 
in Table 5. Flops are complex floating point 
operations[15]. 



Multiplications (flops) 

Autocorrelation 

(2M-p)(l+p)/2 

AR coefficients 
Yule Walker 

P 2 +P 

xpMV coefficients 

P(P+ l)/2 

FFT 

M o.sLog 2 (M os)/ 2 

Total(Approx) 

(M+p)(l+p) 
+M os Log 2 (M os)/ 2 



Additions (flops) 

Autocorrelation 

(2M-p)(l+p)/2 + M 

AR coefficients 
Yule Walker 

X3 

to 

+ 

*3 

ip mv coefficients 

p(p+ l)/2 

FFT 

M osLog 2 (M os) 

Total(Approx) 

(M 4- p)(l 4- p) 4- M 
+M os Log 2 (M os) 


Table 5. Computational effort based on estima- 
tor order p, oversampling os, and the number of 
data points in the range cell M. flops are com- 
plex floating point operations. 

For a typical 2 pm LAWS regime (M = 512, 
Q = 10.0), the maximum optimal order p is 
about 12 and optimal oversampling os is 1. Thus 
the dominant terms for computational effort are 
(2 M - p)(l + p)/2 + MLog 2 (M)/2 for complex 
multiplications and (2 M — p)(l + p)/2 + M + 
M Log 2 (M) for complex additions. 

Computational effort (total of complex mul- 
tiplications and additions) for a typical LAWS 
regime (M=512, Q=7.0) versus estimator per- 
formance is shown in Fig. 2. The solid lines rep- 
resent the total computational effort whereas the 
dashed lines are the computational effort for the 
FFT portion of the estimation routine. For this 
case, computational effort is not dominated by 
the FFT portion of the estimation routine. 

Since computational effort for the CAPON es- 
timator depends on order p, a reduced order re- 
duced performance (reduced from optimal) esti- 
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mator may be of interest for regimes that require 
high optimal estimator order p. A reduced order 
estimator tends to increase the standard devia- 
tion of the good estimates g. 



Figure 2. Estimator performance vs. compu- 
tational effort for M=512 and fi=7.0. Points 
on curve refer to performance for different filter 
order. Solid curves are the total complex oper- 
ations in flops. Dashed curves are the complex 
operations {flops) for the FFT portion of the es- 
timation routine. 

Order reduction is considered for cases where g 
increases no more than 5% and 10%. Perfor- 
mance results due to relaxed estimator order are 
included. 

6 Results M=10 

For this set, performance of the estimator will be 
evaluated for the f2=0.2 and 0.3 cases. For 
larger the signal spectral width becomes com- 
parable to the velocity search space. 

As oversampling os increases, performance ap- 
proaches a constant. Optimal oversampling is 
the minimum oversampling where the perfor- 
mance does not degrade more than the statistical 
variations in the estimates for g. For 12= 0.2 and 
0.3. optimal oversampling os is 3. For the low 
$ cases, oversampling can improve the standard 
deviation of the good estimates g for this pa- 
rameter regime by up to 20%. Performance vs. 
oversampling os for 12=0.2 and 0.3 are shown in 
Fig. 3 and Fig. 4. 

There is a tendency for performance to in- 


crease with order p to the optimal performance, 
then gradually decrease. In addition, optimal or- 
der decreases with increasing <F. Estimator order 
sensitivity for 12=0.2 and 0.3 are shown in Fig. 5 
and Fig. 6. 

Single realizations of the spectral estimates 
can provide insight into the optimal performance 
of an estimator. There is a tendency for the 
width of the signal spectral peak to increase for 
increased 3>. Intuitively one might think that 
a wider spectrum will produce a larger estima- 
tion error g. Results from simulations show that 
at large $ the wider spectral function has less 
error and requires less oversampling for optimal 
performance. Single realizations of the spectral 
estimates for 12=0.2 and 0.3 are shown in Fig. 7 
and Fig. 8. The tick marks represent a sam- 
pling of the spectral estimate with os = 1. The 
realizations demonstrate that the oversampling 
resolves the peak. 

The PDF of the estimators are characterized 
by a clump of good estimates and uniformly dis- 
tributed bad estimates. The PDF’s were mod- 
eled by both a single Gaussian component ( ) 

and a double Gaussian component ( — ) for the 
good estimates. The vertical lines indicate the 
boundary region for estimating the standard de- 
viation of the good estimates g without assuming 
a model for the PDF. The parameters of the PDF 
were estimated using a maximum likelihood algo- 
rithm on the histograms (ML-PDF)[16]. These 
parameters were used to define a region around 
the mean velocity that permits an estimate of g 
without assuming a model for the PDF. 

The bias of the good estimates Av and the 1 a 
error bars are also calculated. For most cases the 
bias is less than 1/1000 of the standard deviation 
of the good estimates g and is not statistically 
different than the error bars. This means there 
is statistically no significant bias for the good es- 
timates. The PDF’s for 17=0.2 and 0.3 are shown 
in Fig. 9 and Fig. 10. Cases of $ = 1, 10, 50, 
and 1000 are reasonable fits with little bias. 

A summary of the performance for the M= 10 
and 12=0.2 and 0.3 is shown in Fig. 11. A sum- 
mary of the optimal order and oversampling is 
shown in Fig. 12. 
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Figure 3. Performance vs. oversampling os for 
M — 10 and fl=0,2. Curves are offset by 0.02. la 
error bars are less than the symbol size. 
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Figure 5. Performance vs. order p for M — 10 
and fi=0.2. Curves are offset by 0.05. la error 
bars are less than the symbol size. 
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Figure 4. Performance vs. oversampling os for 
M — 10 and fl=0.3. Curves are offset by 0.02. la 
error bars are less than the symbol size. 


Figure 6. Performance vs. order p for M — 10 
and fi— 0.3. Curves are offset by 0.05. la error 
bars are less than the symbol size. 4 
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Figure 7. Single realizations of the CAPON Figure 9. PDF of the Velocity Estimates using 
spectral estimates for M = 10 and D=0.2 for the CAPON estimator with M = 10 and D=0.2 
$=1.10.50.1000 with optimal estimator parame- for $=1,10,50,1000. 
t.ers as a function of frequency index k. 





Figure 8. Single realizations of the CAPON Figure i 0 . PDF of the. Velocity Estimates using 
spectral estimates for M = 10 and ft=0.3 for the CAPON estimator with M =10 and «=0.J 
3>= 1,1 0,50, 1000 with optimal estimator parame- f or f or jq 50 !()()() 
ters as a function of frequency index k. 
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Figure 11. Performance of the CAPON estima- 
tor with M= 10 and $7=0.2 and 0.3 using optimal 
estimator parameters, la error bars are less than 
the symbol size. Best fit model is (-). 



Figure 12, Optimal parameters with M~ 10 
and ft=0.2 and 0.3. 


Reducing the order (from optimal) results in 
an increase in the standard deviation g of the 
good estimate. The minimum required estima- 


tor order p with M= 10 and $7=0.2 and 0.3 for no 
greater than 5% and 10% degradation (increase) 
in the standard deviation g of the good estimates 
are shown in Fig. 13. 



Figure 13. Minimum estimator order p with 
M = 10 and J7=0.2 and 0.3 for no greater than 
5% and 10% increase in the standard deviation 
g of the good estimates. 


ft 

%D 

bo 

a 

7 

X 

0.2 

0.00 

7.9381 

.9796 

1.8497 

.7584 

0.3 

0.00 

9.5469 

.9885 

2.2648 

.6478 

0.2 

0.05 

7.8439 

.9846 

1.8447 

.7584 

0.3 

0.05 

9.0815 

.9974 

2.2010 

.6478 

0.2 

0.10 

7.8968 

.9850 

1.8650 

.7584 

0.3 

0.10 

9.9245 

.9827 

2.3456 

.6478 

ft 

%D 

9o 

e 

A 


0.2 

0.00 

199.81 

.6870 

2.3680 

.7554 

0.3 

0.00 

77.824 

.8576 

1.5935 

.6344 

0.2 

0.05 

199.81 

.6870 

2.3680 

.7554 

0.3 

0.05 

77.824 

.8576 

1.5935 

.6344 

0.2 

0.10 

199.81 

.6870 

2.3680 

.7554 

0.3 

0.10 

77.824 

.8576 

1.5935 

.6344 


Table 6. Model Fit Parameters for M= 10, 
$7=0.2 and 0.3. Performance for the optimal and 
no more than 5% and 10% increase in the stan- 
dard deviation g of the good estimates. 
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The performance curves are fitted to empir- 
ical models[10]. The parameters b and g/w v 
are modeled as b($>) = [1 + ( $/6 0 ) a ]~ 7 and 
g/w v {$) = x[l + ($/2o) e ] _< * + M respectively. 
The parameters for optimal (% D=0) and de- 
graded (% D=0.05 and % D=0.10) performance 
are listed in Table 6. Model fits are shown as the 
solid lines for every performance curve. 

7 Results M— 16 

For this set, performance of the estimator will be 
evaluated for the D=0.2, 0.3, 0.5, 1.0 cases. 

Performance vs. oversampling os is shown in 
Appendix A, Fig. A1 - A4. Performance vs. or- 
der p is shown in Appendix A, Fig. A5 - A8. 
Single realizations of the spectral estimates are 
shown in Appendix A, Fig. A9 - A12. The PDF’s 
are shown in Appendix A, Fig. A13 - A16. 

A summary of the performance for the M= 16 
and fl=0.2,0.3,0.5, and 1.0 is shown in Fig. 14. 
A summary of the optimal order and oversam- 
pling is shown in Fig. 15. 



Figure 14. Performance of the CAPON esti- 
mator with M=16 and J7— 0.2, 0.3, 0.5, and 1.0 
using optimal estimator parameters. 1 a error 
bars are less than the symbol size. 



Figure 15. Optimal parameters with M=16 
and fi=0.2, 0.3, 0.5, and 1.0. 



Figure 16. Minimum estimator order p with 
M = 16 and fi=0.2, 0.3, 0.5, and 1.0 for no greater 
than 5% and 10% increase in the standard devi- 
ation g of the good estimates. 
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The minimum required estimator order p with 
M — 16 and fi=0.2 , 0.3, 0.5 and L0 for no greater 
than 5% and 10% degradation (increase) in the 
standard deviation g of the good estimates are 
shown in Fig. 16. 

The model parameters for optimal (% D=0) 
and degraded (% D=0.05 and % D=0.10) per- 
formance are listed in Table 7. Model fits are 
shown as the solid lines in Fig. 14. 
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%D 
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2.5561 

.6362 
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15.860 

1.1228 

4.0381 

.6143 

0.2 

0.05 

7.7098 

1.0484 

1.6617 

.7490 

0.3 

0.05 

8.9319 

1.0710 

2.0077 

.6659 

0.5 

0.05 

11.164 

1.0862 

2.6444 

.6362 

1.0 

0.05 

23.727 

1.0766 

5.5904 

.6143 

0.2 

0.10 

7.6513 

1.0457 

1.6538 

.7490 

0.3 

0.10 

9.3928 

1.0575 

2.0820 

.6659 

0.5 

0.10 

10.879 

1.0885 

2.6020 

.6362 

1.0 

0.10 

24.253 

1.0735 

5.6868 

.6143 

Cl 

%D 

9o 

€ 

6 


0.2 

0.00 

67.5253 

.7032 

1.2938 

.7337 

0.3 

0.00 

59.6000 

.8039 

1.3699 

.6164 

0.5 

0.00 

35.7440 

.8889 

1.3207 

.5004 

1.0 

0.00 

25.8000 

1.1844 

1.2049 

.3832 

0.2 

0.05 

67.5253 

.7032 

1.2938 

.7337 

0.3 

0.05 

59.6000 

.8039 

1.3699 

.6164 

0.5 

0.05 

35.7440 

.8889 

1.3207 

.5004 

1.0 

0.05 

25.8000 

1.1844 

1.2049 

.3832 

0.2 

0.10 

67.5253 

.7032 

1.2938 

.7337 

0.3 

0.10 

59.6000 

.8039 

1.3699 

.6164 

0.5 

0.10 

35.7440 

.8889 

1.3207 

.5004 

1.0 

0.10 

25.8000 

1.1844 

1.2049 

.3832 


Table 7. Model Fit Parameters for M=16, 
f}=0.2, 0.3, 0.5 and 1.0. Performance for the op- 
timal (%D=0.0) and no more than 5% and 10% 
increase in the standard deviation g of the good 
estimates. 

8 Results M=32 

For this set, performance of the estimator will be 
evaluated for the $1=0.2, 0.3, 0.5, 1.0, 2.0, and 
3.0 cases. 


Performance vs. oversampling os is shown in 
Appendix B, Fig. B1 - B5. Performance vs. or- 
der p is shown in Appendix B, Fig. B6 - B10. 
Single realizations of the spectral estimates are 
shown in Appendix B, Fig. Bll - B15. The 
PDF’s are shown in Appendix B, Fig. B16 - B20. 

A summary of the performance for the M —32 
and fi=0.2, 0.3, 0.5, 1.0, 2.0, and 3.0 is shown 
in Fig. 17. A summary of the optimal order and 
oversampling is shown in Fig. 18. 



<D 


Figure 17. Performance of the CAPON esti- 
mator with M = 32 and fl=0.2, 0.3, 0.5, 1.0, 2.0, 
and 3.0 using optimal estimator parameters. 1 o 
error bars are smaller than the symbol size. 

The minimum required estimator order p with 
M = 32 and J7=0.2 , 0.3, 0.5, 1.0, 2.0, and 3.0 for 
no greater than 5% and 10% degradation (in- 
crease) in the standard deviation g of the good 
estimates are shown in Fig. 19. 

The model parameters for optimal (% D=0) 
and degraded (% D=0.05 and % D=0.10) per- 
formance are listed in Table 8. Model fits are 
the solid lines shown in Fig. 17. 
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Figure 18. Optimal parameters with M — 32 
and fi=0.2, 0.3, 0.5, 1.0, 2.0, and 3.0. For Q > 1, 
os — 1 . For Q = 0.2 and 0.5, os— 2. 



Figure 19. Minimum estimator order p with 
M— 32 and 12=0.2, 0.3, 0.5, 1.0, 2.0 and 3.0 for 
no greater than 5% and 10% increase in the stan- 
dard deviation g of the good estimates. 
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1.1434 

1.7283 
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11.832 

1.1541 

2.4022 

.6790 

1.0 

0.05 

18.075 

1.1784 

4.0029 

.6920 

2.0 

0.05 

44.794 

1.1657 

10.196 

.7178 

3.0 

0.05 

127.13 

1.1245 

29.445 

.6506 

0.2 

0.10 

8.2038 

1.1109 

1.5136 

.8411 

0.3 

0.10 

9.0809 

1.1380 

1.7453 

.7203 

0.5 

0.10 

11.940 

1.1494 

2.4213 

.6790 

1.0 

0.10 

18.239 

1.1746 

4.0351 

.6920 

2.0 

0.10 

44.427 

1.1664 

10.121 

.7178 

3.0 

0.10 

127.18 

1.1245 

29.4586 

.6506 

n 

%D 

9o 

e 

8 


0.2 

0.00 

14890.4 

.4027 

9.3678 

.7309 

0.3 

0.00 

84.4809 

.6048 

1.7585 

.6077 

0.5 

0.00 

33.8113 

.7727 

1.3989 

.4927 

1.0 

o 

o 

o 

18.544 

1.0551 

1.0714 

.3664 

2.0 

0.00 

14.1141 

1.3774 

.8169 

.2729 

3.0 

0.00 

18.6496 

1.6891 

.7894 

.2385 

0.2 

0.05 

14890.4 

.4027 

9.3678 

.7309 

0.3 

0.05 

84.4809 

.6048 

1.7585 

.6077 

0.5 

0.05 

33.8113 

.7727 

1.3989 

.4927 

1.0 

0.05 

18.544 

1.0551 

1.0714 

.3664 

2.0 

0.05 

14.1141 

1.3774 

.8169 

.2729 

3.0 

0.05 

18.6496 

1.6891 

.7894 

.2385 

0.2 

0.10 

14890.4 

.4027 

9.3678 

.7309 

0.3 

0.10 

84.4809 

.6048 

1.7585 

.6077 

0.5 

0.10 

33.8113 

.7727 

1.3989 

.4927 

1.0 

0.10 

18.544 

1.0551 

1.0714 

.3664 

2.0 

0.10 

14.1141 

1.3774 

.8169 

.2729 

3.0 

0.10 

18.6496 

1.6891 

.7894 

.2385 


Table 8. Model Fit Parameters for M=32, 
i}=0.2, 0.3, 0.5, 1.0, 2.0, and 3.0. Performance 
for the optimal (%>D=0.0) and no more than 5% 
and 10% increase in the standard deviation g of 
the good estimates. 
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9 Results M=64 

For this set, performance of the estimator will 
be evaluated for the fi=0.3, 0.5, 1.0, 2.0, and 3.0 
cases. 

Performance vs. oversampling os is shown in 
Appendix C, Fig. Cl - C5. Performance vs. or- 
der p is shown in Appendix C, Fig. C6 - CIO. 
Single realizations of the spectral estimates are 
shown in Appendix C, Fig. Cll - C15. The 
PDF's are shown in Appendix C, Fig. C16 - C20. 

A summaxy of the performance for the M — 64 
and ft=0.3, 0.5, 1.0, 2.0, and 3.0 is shown in 
Fig. 20. A summary of the optimal order and 
oversampling is shown in Fig. 21. 

The minimum required estimator order p with 
M = 64 and f)=0.3 , 0.5, 1.0, 2.0, and 3.0 for no 
greater than 5% and 10% degradation (increase) 
in the standard deviation g of the good estimates 
are shown in Fig. 22. 



Figure 20. Performance of the CAPON esti- 
mator with M ~ 64 and f2=0.3, 0.5, 1.0, 2.0, and 
3.0 using optimal estimator parameters. 1 a error 
bars are smaller than the symbol size. 



Figure 21. Optimal parameters with M = 64 
and 17=0.3, 0.5, 1.0, 2.0, and 3.0. For ft > 1, 
os= 1. For 17=0.3 and 0.5, os= 2. 



Figure 22. Minimum estimator order p with 
M — 64 and 17=0.3, 0.5, 1.0, 2.0 and 3.0 for no 
greater than 5% and 10% increase in the stan- 
dard deviation g of the good estimates. 

The model parameters for optimal (% D=0) and 
degraded (% D— 0.05 and % D— 0.10) perfor- 
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mance are listed in Table 9. Model fits are the 
solid lines shown in Fig. 20. 



Table 9. Model Fit Parameters for M=6 4, 
$2 = 0.3. 0.5. 1.0, 2.0. and 3.0. Performance for the 
optimal (no more than 5% and 10% increase in 
the standard deviation g of the good estimates. 

10 Results M— 128 

For this set, performance of the estimator will 
be evaluated for the £7=0.5, 1.0, 2.0, 3.0, and 7.0 
cases. 

Performance vs. oversampling os is shown in 
Appendix D, Fig. D1 - D5. Performance vs. or- 


der p is shown in Appendix D, Fig. D6 - D10. 
Single realizations of the spectral estimates are 
shown in Appendix D, Fig. Dll - D15. The 
PDF’s are shown in Appendix D, Fig. D16 - D20. 

A summary of the performance for the M — 128 
and £7=0.3, 0.5, 1.0, 2.0, and 3.0 is shown in 
Fig. 23. A summary of the optimal order and 
oversampling is shown in Fig. 24. 

The minimum required estimator order p with 
M=128 and $7=0.3 , 0.5, 1.0, 2.0, and 3.0 for no 
greater than 5% and 10% degradation (increase) 
in the standard deviation g of the good estimates 
are shown in Fig. 25. 

The model parameters for optimal (% D=0) 
and degraded (% D=0.05 and % D=0.10) per- 
formance are listed in Table 10. Model fits are 
the solid lines shown in Fig. 23. 



<D 


Figure 23. Performance of the CAPON esti- 
mator with M=128 and $7=0.5, 1.0, 2.0, 3.0, and 
7.0 using optimal estimator parameters. la error 
bars are smaller than the symbol size. 
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Figure 24. Optimal parameters with M = 128 
and Q=0.5, 1.0, 2.0, 3.0, and 7.0. For Q > 1, 
05=1. For £}=0.5, 05=2. 



Figure 25. Minimum estimator order p with 
M = 128 and £2=0.5, 1.0, 2.0, 3.0 and 7.0 for no 
greater than 5% and 10 % increase in the stan- 
dard deviation g of the good estimates. 
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0.5 

0.05 

59.1149 

.6417 

1.8396 

.4927 

1.0 

0.05 

13.2728 

1.2783 

.6908 

.3590 

2.0 

0.05 

10.5849 

1.4301 

.6542 

.2598 

3.0 

0.05 

10.5032 

1.6436 

.5669 

.2142 

7.0 

0.05 

13.4402 

1.7921 

.5521 

.1449 

0.5 

0.10 

59.1149 

.6417 

1.8396 

.4927 

1.0 

0.10 

13.2728 

1.2783 

.6908 

.3590 

2.0 

0.10 
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Table 10. Model Fit Parameters for M= 128, 
J7=0.5, 1.0, 2.0, 3.0, and 7.0. Performance for 
the optimal (%D=0.0) and no more than 5% and 
10% increase in the standard deviation g of the 
good estimates. 

11 Results M=512 

For this set, performance of the estimator will be 
evaluated for the £2=7.0 and 10.0 cases. 

Performance vs. order p is shown in Appendix 
E, Fig. El - E2. Single realizations of the spec- 
tral estimates are shown in Appendix E, Fig. E3 
- E4. The PDF’s are shown in Appendix E, 
Fig. E5 - E6. 
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Figure 26. Performance of the CAPON esti- 
mator with M ~ 512 and 7.0 and 10.0 using 
optimal estimator parameters, la error bars are 
smaller than the symbol size. 



Figure 27. Optimal parameters with A'/=512 
and f]=7.0 and 10.0. For il > 7, os=l. 



Figure 28. Minimum estimator order p with 
A/=512 and Q=7.0 and 10.0 for no greater than 
5 and 10 % increase in the standard deviation of 
the good estimates g. 
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0.05 

14.7845 

2.3502 

.4090 

.1206 

7.0 
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Table 11. Model Fit Parameters for A/= 512, 
12=7.0 and 10.0. Performance for the optimal 
(%D=0.0) and no more than 5% and 10% in- 
crease in the standard deviation g of the good 
estimates. 

A summary of the performance for the A^=512 
and 0=7.0 and 10.0 is shown in Fig. 26. A sum 
mary of the optimal order and oversampling i* 
shown in Fig. 27. 

The minimum required estimator order p with 
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M — 512 and Q=7.0 and 10.0 for no greater than 
5% and 10% degradation (increase) in the stan- 
dard deviation g of the good estimates are shown 
in Fig. 28. 

The model parameters for optimal (% D=0) 
and degraded (% D=0.05 and % D=0.10) per- 
formance are listed in Table 11. Model fits are 
the solid lines shown in Fig. 26. 
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Appendix A, M=16 
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Figure Al. Performance vs. oversampling os 
for M = 16 and 0.2. Curves are offset by 
0.05. (la error bars). 



os 


Figure A3. Performance vs. oversampling os 
for M = 16 and Q— 0.5. Curves are offset by 
0.05. (la error bars). 



Figure A2. Performance vs. oversampling os Figure A4. Performance vs. oversampling os 
for M = 16 and fl=0.3. Curves are offset by for M = 16 and f2=1.0. Curves are offset by 
0.05. (la error bars). 0.05. (la error bars). 
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Figure A5. Performance vs. order p for M = Figure A7. Performance vs. order p for M = 
16 and 11=0.2. Curves are offset by 0.05. 1 a 16 and 11=0.5. Curves are offset by 0.05. lcr 

error bars are smaller than the symbol size. error bars are smaller than the symbol size. 



Figure A6. Performance vs. order p for M = Figure A8. Performance vs. order p for M = 
16 and 11=0.3. Curves are offset by 0.05. lcr 16 and 11=1.0. Curves are offset by 0.05. lcr 
error bars are smaller than the symbol size. error bars are smaller than the symbol size. 
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Figure A9. Single realizations of the CAPON 
spectral estimates for M — 16 and Sl=0.2 for 
$=1,10.50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



Figure All. Single realizations of the CAPON 
spectral estimates for M= 16 and 12=0.5 for 


$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



Figure A10. Single realizations of the CAPON 


spectral estimates for M—16 and 12=0.3 for 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



spectral estimates for M=16 and 12=1.0 fon 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 
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Figure A13. PDF of the Velocity Estimates Figure A15. PDF of the Velocity Estimates 
using the CAPON estimator with M = 16 and using the CAPON estimator with M = lt> and 
0=0.2 for $= 1,10, 50, 1000. 0=0.5 for $=1,10,50,1000. 




Figure A14. PDF of the Velocity Estimates Figure A16. PDF of the Velocity Estimates 
using the CAPON estimator with A/=16 and using the CAPON estimator with M = 16 and 
0=0.3 for for $=1,10,50,1000. 0=1.0 for for $=1,10,50,1000. 
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Appendix B, M=32 



08 


Figure Bl. Performance vs. oversampling os 
for M =32 and 0=0.2. Curves are offset by 0.05. 
(lcr error bars). 



08 


Figure B2. Performance vs. oversampling os 
for M= 32 and 0=0.3. Curves are offset by 0.05. 
(lcr error bars). 



OS 


Figure B3. Performance vs. oversampling os 
for M=32 and 0=0.5. Curves are offset by 0.05. 
(lcr error bars). 



08 


Figure B4. Performance vs. oversampling oa 
for M — 32 and 0=1.0. Curves are offset by 0.05. 
(lcr error bars). 
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08 


Figure B5. Performance vs. oversampling os 
for M = 32 and f2=2.0. Curves are offset by 0.05. 
(la error bars). 



os 


Figure B6. Performance vs. oversampling os 
for M — 32 and f?= 3.0. Curves are offset by 0.05. 
(la error bars). 



P 

Figure B7. Performance vs. order p for M = 32 
and f)=0.2. Curves are offset by 0.05. la error 
bars are less than the symbol size. 
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Figure B8. Performance vs. order p for M = 32 
and f]=0.3. Curves are offset by 0.05. la error 
bars are less than the symbol size. 
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Figure B9. Performance vs. order p for M=32 Figure Bll. Performance vs. order p for 
and (1=0.5. Curves are offset by 0.05. lcr error M= 32 and (1=2.0. Curves are offset by 0.10. 
liars are less than the symbol size. lcr error bars are less than the symbol size. 



Figure BIO. Performance vs. order p for Figure B12. Performance vs. order p for 
M = 32 and (1=1.0. Curves are offset by 0.05. M = 32 and (1=3.0. Curves are offset by 0.10. 
1(7 error bars are less than the symbol size. lcr error bars are less than the symbol size. 
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frequency index k 

Figure B13. Single realizations of the CAPON 
spectral estimates for M — 32 and Q— 0.2 for 
<$>=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



frequency index k 

Figure B14. Single realizations of the CAPON 


spectral estimates for M — 32 and 12=0.3 for 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 
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Figure B15, Single realizations of the CAPON 
spectral estimates for M = 32 and 12=0.5 for 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



Figure B16. Single realizations of the CAPON 
spectral estimates for M = 32 and 12=1.0 fon 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 
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Figure B17. Single realizations of the CAPON 
spectral estimates for M — 32 and 12=2.0 for 
$=1,10. 50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 
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Figure B18. Single realizations of the CAPON 
spectral estimates for M — 32 and 12=3.0 for 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



Figure B19. PDF of the Velocity Estimates 
using the CAPON estimator with M — 32 and 
f2=0.2 for $=1,10,50,1000. 



Figure B20. PDF of the Velocity Estimates 
using the CAPON estimator with M — 32 and 
12=0.3 for for $=1,10,50,1000. 
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Figure B21. PDF of the Velocity Estimates Figure B23. PDF of the Velocity Estimates 
using the CAPON estimator with M = 32 and using the CAPON estimator with M =32 and 
0=0.5 for $=1,10,50.1000. 0=2.0 for $=1,10,50,1000. 



Figure B22. PDF of the Velocity Estimates Figure B24. PDF of the Velocity Estimates 
using the CAPON estimator with M= 32 and using the CAPON estimator with M= 32 and 
0=1.0 for for $=1.10.50,1000. 0=3.0 for for $=1,10,50,1000. 
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Appendix C, M=64 



OS 


Figure Cl. Performance vs. oversampling os 
for M— 64 and ft— 0.3. Curves are offset by 0.05. 
1 a error bars are smaller than the symbol size. 



os 


Figure C3. Performance vs. oversampling os 
for Af= 64 and fl=1.0. Curves are offset by 0.05. 
la error bars are smaller than the symbol size. 



Figure C2. Performance vs. oversampling os Figure C4. Performance vs. oversampling os 
for M=64 and 17=0.5. Curves are offset by 0.05. f or M= 64 and Vt= 2.0. Curves are offset by 0.05. 
la error bars are smaller than the symbol size. ( \ G error bars). * 
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08 


Figure C5. Performance vs. oversampling os 
for M =64 and $1=3.0. Curves are offset by 0.05. 
(la error bars). 



Figure C6. Performance vs. order p for M= 64 
and $1=0.3. Curves are offset by 0.05. la error 
bars are less than the symbol size. 



Figure C7. Performance vs. order p for M ~ 64 
and D=0.5. Curves are offset by 0.05. la error 
bars are less than the symbol size. 



Figure C8. Performance vs. order p for M= 64 
and $1=1.0. Curves are offset by 0.05. la error 
bars are less than the symbol size. 
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Figure C9. Performance vs. order p for M = 64 
and $7=2.0. Curves are offset by 0.05. la error 
bars are less than the symbol size. 



Figure CIO. Performance vs. order p for 
M ~ 64 and $7=3.0. Curves are offset by 0.10. 
la error bars are less than the symbol size. 



frequency index k 

Figure Cll. Single realizations of the CAPON 
spectral estimates for Af= 64 and $7=0.3 for 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



frequency index k 

Figure C12. Single realizations of the CAPON 
spectral estimates for M=64 and $7=0.5 foi 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 
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Figure C13. Single realizations of the CAPON 
spectral estimates for M — 64 and 52=1.0 for 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



frequency index k 

Figure C14. Single realizations of the CAPON 
spectral estimates for M = 64 and 12=2.0 for 
$ = 1,10.50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 
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Figure C15. Single realizations of the CAPON 
spectral estimates for M = 64 and 12=3.0 for 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k . 



Figure C16. PDF of the Velocity Estimates 
using the CAPON estimator with M =04 and 
12=0.3 for for $=1,10,50,1000. 
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Figure Cl 7. PDF of the Velocity Estimates Figure C19. PDF of the Velocity Estimates 
using the CAPON estimator with M = 64 and using the CAPON estimator with M = 64 and 
D=0.5 for $=1.10,50,1000. fl=2.0 for $=1,10,50,1000. 




Figure C18. PDF of the Velocity Estimates Figure C20. PDF of the Velocity Estimates 
using the CAPON estimator with M=64 and using the CAPON estimator with A/=64 and 
D=1.0 for for $ = 1,10,50,1000. fi= 3.0 for for $=1,10,50,1000. 
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Appendix D, M=128 




Figure Dl. Performance vs. oversampling os Figure D3. Performance vs. oversampling os 
for M = 128 and fi=0.5. Curves are offset by 0.05. for M = 128 and 0=2.0. Curves are offset by 0.05. 
(lcr error bars). (l<r error bars). 
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Figure D2. Performance vs. oversampling os 
for Af =128 and 1.0. Curves are offset by 0.05. 

Figure D4. Performance vs. oversampling os 
for Af= 128 and Q=3.0. Curves are offset by 0.05. 

(lcr error bars). 

(la error bars). 











Figure D5. Performance vs. oversampling os Figure D7. Performance vs. order p for 
for M=128 and f}=7.0. Curves are offset by 0.05. M= 128 and fi— 1.0. Curves are offset by 0.10. 
lcr error bars are less than the symbol size. lcr error bars are less than the symbol size. 



Figure D6. Performance vs. order p for Figure D8. Performance vs. order p for 
M=128 and f2=0.5. Curves are offset by 0.10. M = 128 and f)=2.0. Curves are offset by 0.10. 
la error bars are less than the symbol size. la error bars are less than the symbol size. 
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Figure D9. Performance vs. order p for 
M — 128 and $1=3.0. Curves are offset by 0.10. 
la error bars are less than the symbol size. 



Figure DIO. Performance vs. order p for 
M — 128 and 12=7.0. Curves are offset by 0.10. 
la error bars are less than the symbol size. 



Figure Dll. Single realizations of the CAPON 
spectral estimates for A/=128 and 12=0.5 for 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



frequency index k 

Figure D12. Single realizations of the CAPON 
spectral estimates for M=128 and $1=1.0 fof 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 
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Figure D13. Single realizations of the CAPON 


spectral estimates for M~ 128 and ft=2.0 for 
$—1,10.50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



Figure D15. Single realizations of the CAPON 
spectral estimates for M — 128 and ft=7.0 for 
$=1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



Figure D14. Single realizations of the CAPON 
spectral estimates for M — 128 and ft =3.0 for 
$ = 1,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 


| * T ■ r — 

: #=10 

— I 1 — 7 — 1 — ; 

0=0.34083 ; 

-0*0.5 

4-0.00143 

' M=128 

9 fw =0.854 

ord=42 

<-0.003 

M =2 

Avrt* =0.007 


4-0.004 


9 % *0.841 


(ML-POF) : 

— l 4-i 1 


-t ■ t ' I t 
0=1 

1 1 1 • 
b=0. 54307 ■ 

0*0.5 

<-0.00081 

M=128 

*1.034 ' 

ord=71 

4-0.012 

o «=2 

dv/w *0.026 . 


+-0.013 


g/w = 1.100 


(ML-POF) 

« * 

' 1 * 1 

« • 


i * i * ■ u ■ ■ ■ ■ i 

0.0 0.2 0.4 0.6 O.t to 


■0=1000 

r M=128 

ord=6 

ot *2 

r 

r 

I i | » — j 

bxO.OOOOO : 
4 - 0.00000 
g/w *0.510 i 
4-0.001 

Av/w =0.000 ’ 
4-0.002 ] 

] 

g/w *0.510 
(ML-POF) ^ 

] 

0=50 

br0.02156 : 

U=0.5 

4*0.00049 

. M=126 

gto =0.684 „ 

:ocd*27 


01=2 

Av/W =0.003 : 


<-0.002 

1 


, gft* =0.683 1 



• 6 

* * ' 1 1 

i i i i 


J l—6—l k i 111 I 1 1 1 

10 0.0 0.2 0.4 0.6 0.1 to 


v/v 


Figure D16. PDF of the Velocity Estimates 
using the CAPON estimator with M = 128 and 1 
0=0.5 for for 0=1,10,50,1000. 
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Figure D17. PDF of the Velocity Estimates Figure D19. PDF of the Velocity Estimates 
using the CAPON estimator with M = 128 and using the CAPON estimator with M= 128 and 
0=1.0 for $=1,10,50,1000. 0=3.0 for $=1,10,50,1000. 




Figure D18. PDF of the Velocity Estimates Figure D20. PDF of the Velocity Estimates 
using the CAPON estimator with M = 128 and using the CAPON estimator with M=128 and 
0=2.0 for for $=1,10,50,1000. 0=7.0 for for $=1,10,50.1000. 
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Appendix E, M=512 



p 


Figure El. Performance vs. order p for M = 512 
and f2=7.0. Curves are offset by 0.05. (lcr error 
bars). 



P 


Figure E2. Performance vs. order p for M=512 
and ft=10.0. Curves are offset by 0.05. (la error 
bars). 



frequency index k 

Figure E3. Single realizations of the CAPON 
spectral estimates for M — 512 and Q=7.0 for 
$=2,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 



frequency index k 

Figure E4. Single realizations of the CAPON 
spectral estimates for M=512 and fl=10.0 foi* 
$=2,10,50,1000 with optimal estimator parame- 
ters as a function of frequency index k. 
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Figure E5. PDF of the Velocity Estimates 
using the CAPON estimator with M — 512 and 
Q=7.0 for $=2,10,50,1000. 
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Figure E6. PDF of the Velocity Estimates 
using the CAPON estimator with M= 512 and 
0=10.0 for for $=2,10,50.1000. 
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